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Abstract 



A new possibility for solving the strong CP-problem is suggested. It is based on the 
concept of a mirror world of particles, with the gauge symmetry and Lagrangian completely 
identical to that of the observable particles. We assume that the ordinary and mirror sectors 
share the same Peccei-Quinn symmetry realized a la Weinberg- Wilczek, so that the #-terms 
are simultaneously canceled by the axion VEV in both worlds. This property remains valid 
even if the symmetry between two sectors is spontaneously broken and the the weak scale 
of the mirror world is larger than the ordinary weak scale, in which case also the mirror 
QCD scale becomes larger than the ordinary one. In this situation our axion essentially 
represents a Weinberg- Wilczek axion of the mirror world with quite a large mass, while it 
j^ ' couples the ordinary particles like an invisible axion. The experimental and astrophysical 

limits are discussed and an allowed parameter window is found with the Peccei-Quinn scale 
f a ~ 10 4 — 10 5 GeV and the axion mass m a ~ 1 MeV, which can be accessible for future 
experiments. We also show that our solution to the strong CP-problem is stable against 
the Planck scale induced effects. 
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1 Introduction 

The strong CP problem remains one of the puzzling points of the modern particle physics (for 
a review, see e.g. ref. [|J). It is related to the P- and CP- violating term 9GG in the QCD 
Lagrangian which is contributed also by the complex phases in the quark mass matrices Mu,d > 
so that its effective value becomes 9 = 9 + arg(det My det Mp). The CP-violating phenomena 
in weak interactions indicate that the CP-violating phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix is large, and hence the quark mass matrices should contain large complex 
phases. Therefore, one would generally expect that 9 is order 1. However, the experimental 
bound on the neutron electric dipole moment yields 9 < 10~ 10 . The strong CP problem simply 
questions why the effective #-term is so small, or in other words, what is the origin of such a 
strong fine tuning between the initial value 9 and the phases of the quark mass matrices? 

The most appealing solution to the problem is provided by the Peccei-Quinn (PQ) mecha- 
nism [|j, based on the concept of spontaneously broken global axial symmetry C/(1)pq. In the 
context of this mechanism 9 essentially becomes a dynamical field, 9 = a/f a , with an effective 
potential induced by the non-perturbative QCD effects which fixes vacuum expectation value 
(VEV) at zero. Here a stands for the pseudo-Goldstone mode of the spontaneously broken 
PQ symmetry, an axion, and f a is a VEV of scalar (or a VEV combination of several scalars) 
responsible for the C/(1)pq symmetry breaking. This scale is model dependent. In particular, 
in the original Weinberg-Wilczek (WW) model |3J f a is order electroweak scale, but in the 
invisible axion models |I], [5], |(| it should be much larger. 

In the most general context, the axion couplings to fermions and photons are described by 
the following Lagrangian terms: 

£a = iQaka^kl^k + -^-aF^F^ + .... (1) 

where the coupling constants are essentially determined by the PQ scale f a . Typically, up to 
model dependent order one coefficients, we have g a k ~ fa^ m ki where m^ is the fermion mass 
(e.g. TTifc = m e ,mN,--- for electron, nucleons, etc.), and g ai ~ (a/2ir)f~ 1 , with a being the 
fine structure constant. 

As for the axion mass, it also depends on the QCD scale A where strong interactions 
become nonperturbative. In the theory with no light quarks (with masses below A), the axion 
mass would emerge due to instanton-induced potential which can be computed in the dilute 
gas approximation, P{a) m —Kcos(a/f a ) with K ~ A 4 , and so we would have m 2 a ~ A 4 //^. 
However, in reality the axion state a gets small mixing with the pseudo-Goldstone bosons it, r], 
etc., related to the chiral symmetry breaking by the light quark condensates V = (qq) ~ A 3 . 
In order to find the true mass eigenvalue corresponding to the axion, the total mass matrix 
should be diagonalized. In doing so, one finds a general expression J?]]: 

m 2 _ El V J£_ (2) 

a fl V + KTtM- 1 { ' 

where N stands for the color anomaly of C/(1)pq current,^] and M = di&g(m u ,md, ■■■) is a mass 
matrix of light quarks, with m q < A. Thus, the axion mass is given as rr? a ~ m q A 3 / f%. More 
precisely, neglecting the s-quark contribution and using the relation (m u + m^iqq) = rn^f 2 , 

4 The PQ charges are normalized so that each of the standard fermion families contributes as N — 1. Hence, 
in the WW model we have iV = N g , where N g (= 3) is the number of fermion families. The same holds in the 



Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) model M. Other models of the invisible axion, e.g. the hadronic 
axion H or archion [H], generally contain some exotic fermions and so iV 7^ N g . 



from (g) one directly arrives to the more familiar formula: 

N(m u m d V\ 1 ' 2 Nz 1 / 2 Um n /lO 6 GeV\ Ar 

m a = -7- = • — : — ~ — - jt-l x 6.2 ev (3) 

fa\m u + mj 1 + Z f a \ f a /N ) K > 

where z = m u /md ^ 0.57, and m,, f n respectively are the pion mass and decay constant. 

In the WW model || the PQ symmetry is broken by two Higgs doublets H\p with the 
VEVs t>i t 2, namely f a = (v/2) sin 2/3, where u = (uf 4- t"!) 1 ^ 2 — 247 GeV is the electroweak 
scale and tan/3 = 1)2/ V\. Therefore, the WW axion is quite heavy, and its mass can vary from 
few hundred keV to several MeVs: 

ww 2N /m u l/y/ 2 150keV 



vsin2/3 U + W sin2/3 ' V ; 

However, its couplings to fermions and photon are too strong and by this reason the WW 
model is completely ruled out for any values of the parameter j3 by a variety of the terrestrial 
experiments as are the search of the decay K + — ► ir + a, J/ip and T decays into a + 7, nuclear 
deexcitations via axion emission, the reactor and beam dump experiments, etc. ||. 

For a generic invisible axion, these experimental data imply a lower bound f a > 10 4 GeV 
or so, which leads to the upper limit on m a of about 1 keV. However, for an axion being so 
light, the stringent bounds emerge from astrophysics and cosmology. In particular, for the 
DFSZ model the combined limits from the stellar evolution and the supernova neutrino signal 
exclude all scales f a up to 10 10 GeV jl], [| , which in turn implies m a < 10~ 4 eV.(] On the other 
hand, the cosmological limits related to the primordial oscillations of the axion field or to the 
non-thermal axion production by cosmic strings, demand the upper bound f a < 10 10-11 GeV 
[p]. |J. Thus, not much parameter space remains available. 

In addition, the invisible axion models have a naturality problem related to the Planck 
scale induced terms. Many arguments suggest that the non-perturbative quantum gravity 
effects do not respect the global symmetries |l(J and thus they can induce the higher order 
effective operators cutoff by Mp which explicitly violate the PQ symmetry. For large f a these 
explicit terms can be very dangerous for the stability of the PQ mechanism. For example, in 
the DFSZ model where L^(1)pq is broken by a gauge singlet scalar 5, the operator -^ — \- h.c. 



would induce 8 > 10 10 unless f a < 10 GeV [10|. Therefore, the Planck scale induced effects 



leave no room for the invisible axion models. On the other hand, the WW model does not 
have this problem and neither it is sensitive to astrophysical constraints, but it is excluded by 
the laboratory limits because of too strong couplings with the matter |j 

In the present paper we suggest a new model for the axion. We assume that there exists 
a parallel sector of "mirror" particles [jl3| , [14| , |15| , 16 1 which is completely analogous to the 



ordinary particle sector, i.e. it has the same gauge group and all coupling constants are equal 
in two sectors. In other words, the whole Lagrangian is invariant with respect a discrete mirror 
parity (M-parity) under the exchange of two sectors. We further assume that the ordinary and 
mirror worlds have the common PQ symmetry realized d la Weinberg- Wilczek, with the £/(1)pq 
charges carried by the ordinary Higgs doublets H\^ and their mirror partners H[ 2 \\ 

5 In the case of the hadronic axion H or archion |6[, a small window around f a ~ 10 6 GeV and so the axion 
mass of few eV can be also permitted. 

6 There are also "axionless" models where the solution to the strong CP-problem is due to spontaneously 
broken P- or CP-invariance fill]. These scenarios, however, have many specifics and do not work in the general 
context. The Planck scale induced effects for the axionless models have been studied in ref. 113] . 

7 The possibility to use mirror symmetry for solving the strong CP problem, in different context, was first 
considered by Rubakov |L7[ . We will briefly comment on this model later. In principle, two sectors could have 
also other common symmetries, e.g. the flavour symmetry |B3. 



If M-parity is an exact symmetry, then the particle physics should be exactly the same in 
two worlds: in particular, the initial #-terms are equal, 6 = 6', the quark mass matrices are 
identical, Mu,d = M{j D , the QCD scales coincide, A = A', and the axion couples both QCD 
sectors in the same way: f~ 1 a(GG + G'G'), so that their non-perturbative dynamics should 
produce the same contributions to the axion effective potential. Clearly, in this case the strong 
CP-problem is simultaneously solved in both worlds - the axion VEV cancels the #-terms both 
in the mirror and ordinary sectors. In such a realization, however, f a remains order 100 GeV 
and thus it is excluded by the same phenomenological grounds as the original WW model. 

As it was suggested in ref. jjlj|, the M-parity can be spontaneously broken and the elec- 
troweak symmetry breaking scale v' in the mirror sector can be larger than the ordinary one 
v. This would lead to somewhat different particle physics in the mirror sector and it is not 
a priori clear that the strong CP-problem still can be simultaneously fixed in both sectors. 
However, we show below that this is just the case! The simple reason is that the softly broken 
M-parity leads to the VEV difference between the ordinary and mirror Higgs doublets, but the 
structure of the Yukawa couplings remains the same in two sectors. 

As far as U(1)pq is a common symmetry for two sectors, now the PQ scale is determined 
by the larger VEV v' , i.e. f a ~ (v f /2) sin 2/3'. The axion state a dominanly comes from the 
mirror Higgs doublets H[ 2 , up to small (~ v/v') admixtures from the ordinary Higgses H±2- 
Hence, it is a WW-like axion with respect to mirror sector, while it couples the ordinary matter 
as an invisible DFSZ-like axion. Therefore, the experimental limits on the axion search require 
f a > 10 4 GeV or so. 

However, there are good news concernin the axion mass - for a given scale f a it appears 
to be much larger than the mass of the ordinary DFSZ axion (|3|). The reason can be briefly 
explained as follows. Due to the larger weak scale, v' ^> v, the mirror quark masses are scaled 
up by a factor v'/v with respect to the ordinary quarks, m' u d ~ (v' /v)m u> d,.„- By this reason, 



also the mirror QCD scale becomes somewhat larger than the ordinary one: A' > A [15, 16 1. 
Therefore, the axion potential dominantly emerges from the mirror non-perturbative dynamics 
and we have m a ~ m! A' 3 //^ or m^ ~ A' 4 //^ (depending on whether the mass of lightest mirror 
quark m' is smaller or larger than A'). Therefore, our axion can be much heavier than the 
ordinary invisible axion with m\ ~ m, g A 3 //^. In particular, if with increasing v' also A' grows 
fastly enough, our axion can maintain the weight category of the WW axion (||) and its mass 
can remain in the MeV range even for f a > 10 4 GeV. 

The paper is organized as follows. In sect. 2 we present our model and show how it solves 
the strong CP-problem. We also compute the axion mass and its couplings to the fermions and 
photons. In sect. 3 we confront our model with the experimental and astrophysical bounds 
and demonstrate that there is an allowed parameter space which can be of interest for the 
future experimental search. The stability of our model against the Planck scale induced effects 
is discussed in sect. 4. Finally, in sect. 5, we summarize our results. 

2 Mirror World and the Peccei-Quinn symmetry 

It was suggested a long time ago that there can exist a parallel (mirror) gauge sector of particles 
and interactions which is the exact copy of the visible one and it communicates the latter 
through the gravity and perhaps also via some other interactions. Various phenomenological 
and astrophysical implications of this idea have been studied in refs. ||13| , |l4| , 15, 16 1. 



In particular, one can consider a model based on the gauge symmetry G x G' where G = 
SU(3) x SU(2) x U{\) stands for the standard model of the ordinary particles: the quark 
and lepton fields ifii = qi, k, uf, d^, ef (i is a family index) and two Higgs doublets Hip, 



while G' = SU(3)' x SU(2)' x U(l)' is its mirror gauge counterpart with the analogous particle 
content: the fermions tp^ = q' { , l' { , uf, df, ef and the Higgses H[ 2 . From now on all fields 
and quantities of the mirror sector have an apex to distinguish from the ones belonging to the 
ordinary one. All fermion fields ijj, ip' are taken in a left-chiral basis. 

Let us assume that the theory is invariant under the mirror parity M: G «-> G' , which 
interchanges all corresponding representations of G and G' . Therefore, the two sectors are 
described by identical Lagrangians and all coupling constants (gauge, Yukawa, Higgs) have the 
same pattern in both of them. In particular, for the Yukawa couplings 

£ Yuk = G%u\q d H 2 + G%d\ qj B x + G%e\l 3 H x + h.c, 

£ Yuk = G'^^H, + G'ijdfq'^ + G'Hefl'jHi + h.c. (5) 

we have Gy D E = Gff D E . In addition, also the initial #-terms are equal, 6 = 9'. 

We further assume that two sectors have a common Peccei-Quinn symmetry U(1)pq under 
which fermions i/ji^'i change their phases by a factor exp(— iui/2) and the Higgses Hi j2 ,H[ 2 
change phases as exp(zo>). Then the renormalizable Higgs Potential has a general form Vtot = 
V + V + V m ix, where for the ordinary Higgses we have: 

V = -hIh{h x - ii\H\H 2 + AiOffJffi) 2 + \2(H\H 2 f + A(if 1 i/ 2 )t(if 1 if 2 ), (6) 

the mirror Higgs potential V has exactly the same pattern by M-parity: H\ 2 — » H[ 2 , and the 
mixing terms are: 

V mix = -KiHiHzfiHiHb + h.c. (7) 

where the coupling constant n should be real due to M-parity|] 

The mirror parity can be spontaneously broken so that the weak interaction scales are 



different in two sectors [15, 16]. The easiest way is to introduce a real singlet scalar rj which is 



odd under M-parity: r\ — ► — rj [18] . Then the following renormalizable terms can be added to 
the Higgs potential: 

AV = h(r ] 2 -fi 2 ) 2 + a k r ] \HlH k + H^H' k ) + p k i2r ] (HlH k -H^H' k ), k = 1,2 (8) 

A non-zero VEV (rj) = jjl spontaneously breaks the M-parity and induces the different effective 
mass 2 terms for the ordinary and mirror Higgses: —^\i e ^\ = — n\ + <?k^ 2 ' + Pk^ 2 and — /J> k 2 ( e g\ = 
—jjl\ + akjJ. 2 — PkP 2 - Therefore, the mirror VEVs v[ 2 are different from the ordinary ones v\ j2 , 
and in general also the VEV ratios v 2 /v\ = tan/3 = x and v' 2 /v[ = tan/3' = x' must be different 
in two sectors. In particular, we will be interested in a situation when v [ 2 3> vi 2 , i.e. the 
weak interaction scale v' = (v 12 + v^ 2 ) 1 ' 2 in the mirror sector is much larger than the ordinary 
weak scale v = (v 2 + v^) 1 / 2 = 247 GeV, v'/v = ( 3> 1 [|l|, |l6|. Clearly, the large hierarchy 



between v' and v requires the fine tuning of the parameters in the Higgs potential. We do not 



discuss here this question and refer the reader e.g. to ref. [16]. 

As far as the ordinary and mirror Yukawa constants in (|5|) are the same, the quark and 
lepton mass matrices essentially have the same form in both sectors: Mjj = Gu{H 2 ), M\j = 
Gu(H 2 ), Md = Gd{H\), M' d = Gd{H[), etc. Thus, the fermion mass and mixing structure 
in the mirror sector should be the same as in the ordinary one, with the mirror up quark 
masses scaled as m' u ct = C,2mu,c,t and the down quark and charged lepton masses scaled as 

8 Notice, that in the limit k = there emerge two separate axial global symmetries, U(1)a for ordinary sector 
under which ipi — > exp(— iw/2)ipi and Hi t 2 —^ exp(io;)i/i i 2, and U(1)'a for mirror sector: tp'j — > exp(— iu}'/2)lp'i 
and H[ 2 ~ * exp(ia»')-ff li2 - Therefore, the term Vmix demands that lo' = lo and thus it reduces U(1)a x U(1)' a 
to its diagonal subgroup U(1)pq. 
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Figure 1: Renormalization group evolution of the strong coupling constants a s and a' s from higher to 
lower energies, for v' = 10 6 GeV and x = x' . 



m 'dsb = d m d,s,b, m 'e,v,T = Cim e ,ix,r, where C2 = v' 2 /v 2 = C(sin/3'/sin/3) and Ci = v'Jvx = 
C(cos/?7cos/3). 

At very high energies, \i 3> v', the strong coupling constants a s (fx) and a' s (fi) should be 
equal due to M-parity. Under the renormalization group (RG) evolution they both evolve down 
in parallel ways until the energy reaches the value of mirror-top mass m' t ~ Qi'mt- Below it 
a' s will have a different slope than a s , and this slope will change every time below the mirror 
quark thresholds m' b ~ Ci m f>> etc, In the evolution of a s these thresholds occur at lower scales, 
fi = mt,mb, etc. Then it is very easy to determine the scale A' at which a' s becomes large, 
once we know that for the ordinary QCD this happens at A ~ 200 MeV. In other words, A' 
becomes a function of (1 2, and for v' ^> v one could obtain a significant difference between the 
QCD scales: A' > A. 

Let us consider for simplicity the case x = x', which yields Ci = C2 = v'/v. Then for the 
values of v' up to 10 4 GeV or so one obtains an approximate scaling low: 



A 



i>2-i>e 



0.28 



0) 



where b n = 11 — ^n are the RG coefficients for the case of n light quarks. For v' > 10 4 
GeV the mirror light quarks masses become larger than A', since they grow faster with v': 



in, 



i/ m u,d = v'/v, and the approximation (||) is not valid - for fi < 
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u,d 



the mirror QCD 



becomes the pure gluodynamics and hence A'/A starts to increase faster with v' (see Fig. 1). 
Let us discuss now the axion physics in our model. Since we are interested only in Goldstone 
modes, let us express the neutral components of the ordinary and mirror Higgs fields as: 



Hi = -y=exp(i(f> k /v k ), 



H 



/0 



-sfflPW* 



V2 



1,2 



(10) 



i>l : 2 and 
-1/ 



1 2 which corresponds to the axion field, 
and ^ - -'- 1 ' 



and identify the combination of the phases 

We immediately observe that the states <pz = v~ l (v2(p2 — ^l^i) and <p' z = V~'{v' 2 <p2 ~ ^iv^i 
are the true Goldstone modes eaten up respectively by the standard and mirror gauge bosons 
Z and Z' which become massive particles by the Higgs mechanism. 

In the limit k = 0, two remaining combinations <f> = v~ 1 (vi4>2 + ^20i) and </>' = f /_1 (u^2 + 
v 2 (f>i) would represent the Goldstone modes of larger global symmetry U(1)a x U(1)' a , respec- 
tively with the decay constants / = vyv^jv = (v/2)sin2/3 and /' = v[v 2 /v' = (v'/2) sin 2/3'. 



Thus, we would have two axion states separately for the ordinary and mirror sectors, and their 
masses would emerge respectively due to quantum anomalies of the U(1)a and U{\)' A currents 
related to ordinary and mirror QCD. In particular, <p would represent the familiar WW axion 
with order ~ 1 MeV mass given by eq. (Q). 

However, the mixing term (Q) reduces the larger global symmetry U(1)a x U(V)' a to a 
common PQ symmetry £7(1)pq. Substituting ( ^0|) in V m j x , we obtain: 



, , ( (p <P \ KVV'.,, ,. , 2 



V mix = -^vWA cos [j-jj)= ffjjif'^ - m 2 + ■■■ (11) 

Thus, the combination A = f~ l (f'cf) — f (/)'), where f a = (f 2 + f' 2 ) 1 ' 2 , receives a mass 2 term: 

which is positive for k > 0. The present experimental limits on the Higgs search imply the 
lower bound on Ma about 100 GeV ||. As for the other combination a = fa l (f4> + /'</>'), it 
remains as a Goldstone mode of the C/(1)pq symmetry to be identified with the axion field. 

Let us consider now the effective #-terms in two sectors. Recalling that 6 = 6' and Gjj,d = 
G'u D due to M-parity, we obtain:^ 

6 = 6 + arg(det M v det M D ) =6 + arg(det G v det G D ) + N(<p/f) 

6' = 6' + arg(det M' v det M' D ) = 6 + arg(det G v det G D ) + N (</>'//') (13) 

where the last terms describe contributions of the phases in ((H]), and they can be immediately 
rewritten as </>'/f = f- x [a-{f / f)A] and 0// = Z" 1 [a +(/'//) A]. We see that the axion field 
a contributes both eqs. ( |l3| ) in the same way, whereas the heavy state A is irrelevant since it 
has a vanishing VEV. Therefore, 6 and 6' will be simoultaneously cancelled by the axion VEV 
fa 1 {a) and so in our model the strong CP-problem is solved in both sectors j^j 

In addition, for the case /' S> / we obtain a quite interesting situation. The axion state a 
dominantly emerges from the phases of the mirror Higgses H[ 2 while it contains only a small 
admixture from the ordinary Higgses H\^'- a m </>' — (f / f ')(/). Therefore, it appears to be a 
WW-like axion with respect to the mirror sector, while from the view of the ordinary sector it 
is just a DFSZ-like invisible axion with a decay constant f a ~ /' = (v'/2) sin 2 /3'.P| 

However, there is a dramatic difference from the DFSZ case. Since the mass of our axion 
receives contributions from both ordinary and mirror QCD, in a full analogy to eq. (||) one 
can write: 

m2 = "1 ( YE + Y^L ^ (14) 

fa \V + KTiM- 1 V' + K'TvM'- 1 J { ' 

where M = diag(?7i( t , m' d ) is the mass matrix of the mirror light quarks, and the values K' ~ A' 4 
and V ~ A' 3 characterize the mirror gluon and quark condensates, with A' being the QCD scale 
in the mirror sector. The first term in this expression gives an ordinary small contribution (|3|) 
and it is negligible in comparison to the second term which can give much larger contribution 

9 In the following, we maintain parametrically the number of fermion families, N = N g , in general formulas, 
but for the numerical calculations we always take iV = 3. 

10 Actually 8 can get radiative corrections due to the fact that the electroweak scales are different in two 
sectors. However, it is well-known that these corrections are negligibly small [h9[. 

11 This is very easy to understand from the following observation. The combination H[H'2 = S in W\) is a 
standard model singlet carrying the double PQ charge, and it can be considered as a composite operator which 
communicates the standard Higgses Hi t 2 via the term V m i x = (H1H2) S + h.c, just as in the DFSZ model. 



as far as m' u d 3> m U) d and A' > A. Certainly, the exact form of the latter term depends on how 
many light quarks (with masses less than A') are contained in mirror sector. In particular, for 
the cases of two (v! and d'), one («'), or zero light quarks, the second term in brackets of eq. ( |l4|) 
becomes respectively (1 + z')^ 1 ^^', m' u V or K', where z' = m' u /m' d = (C2/Ci) z = (x'/x)z. 

Let us take for example the case x = x', and consider the regime when at least u' is light, 
m' u = (v'/v)m u < A' (recall that this regime holds up to the values v' ~ 10 4 GeV or so). Then 
for the axion mass holds an approximate scaling low: 

m ^V^w{TT^) "{a) WJ ma ■ (15) 

with respect the ordinary WW axion mass (0). Using also eq. (||), we thus obtain: 

fv\ om 150 keV , , 

m a - -7 x . na , (16) 



sin 2/3 

In general case, for very large v' and arbitrary x', the simple scaling low ( |l5|) is no more valid, 
since for x ^ x' the masses of u' and d' scale differently and in addition, with increasing v', 
only one or none of them may remain light, with mass less than A'). In this situation one has 
to use the general formula fll4|). Nevertheless, m a still decreases very slowly with increasing v' 
and for enough large x, x', its value can be of the order of MeV. 

The results of numerical computations from the general formula ( |l4| ) for some interesting 
cases are presented in Fig. 2. The relevant parameters are taken as x = tan/3, x' = tan/3' 
in the range 1 — 100, and f a = v'/(x' + x'^ 1 ). We use one-loop RG evolution for the strong 
coupling constants with A = 200 MeV. The following values of the quark masses are used: 
m u = A MeV, m d = 7 MeV, m s = 150 MeV (at /jl = 1 Gev), and m c (m c ) = 1.3 GeV, m b = 4.3 
GeV and mt = 170 GeV (respectively at /x = m c ,mb,mt)- For the parameters V and K 
related to the quark and gluon condensates we take V = (250 MeV) 3 and K = (230 MeV) , 
and assume that the similar parameters in the mirror sector scale as V'/V = (A'/A) 3 and 
K'/K = (A'/A) 4 . We see that for enough large values of x,x' the axion mass can be order 1 
MeV even for /„ = 10 4 - 10 6 GeV. 

Let us discuss now the interaction terms (|l|). Clearly, our axion couples the ordinary 
fermions and photons in exactly the same way as the DFSZ axion with the PQ scale f a . 
Hence, for the Yukawa coupling constants to the electron, down quark and up quark we can 
use the familiar expressions |lj: 



— sin /3, g ad = — 

Ja Ja 



■ 2R N 

sin p 



l + z 
while for the axion-photon coupling we have: 



_ m u 

9au — „ 
Ja 



cos 2 /3 



Nz 

l + z 



(17) 



^^iVC _ 8 6KTr(M~ 1 Q 2 ) _2z_ 

9ai ~ 2vr f a ' 3 V + KTr(Af-i) " 1 + z [ j 

where the trace is taken over the light quark states (u, d), Q are their electric charges (+4/3, —1/3), 
and a = 1/137 is the fine structure constant. In addition, our axion couples to mirror photons, 
with the constant: 

' ^^^ _ 8 6K'Tr(M'- l Q' 2 ) 

9a ^ 2vr f a ' 3 V + K'Tt(M'- 1 ) l ' 

where the factor C for the case of two (u',d'), one (u') or no light quarks respectively takes 
the values 2z'/(l + z'), and 8/3. 
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Figure 2: Contour plots for the axion mass (solid) as a function of f a and x' , for the cases x = x' 
and x = 60. The regions left from the dash-dots curves are excluded by the observed neutrino signal 
from SN 1987A while the regions right from the dash curves are excluded by the SMM bounds on the 
gamma ray flux from SN 1987A. The vertical dott line indicates a conservative bound on the PQ scale 
from the terrestrial experiments. 



Hence, the axion decay widths into the visible and mirror photons respectively are: 



T(a -> 77) 
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In addition, if m a > 2m e , the axion can decay also into electron-positron pair: 



r(a->e + e _ ) 



9le m a 
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3 Experimental and astrophysical bounds 

Let us discuss now the constraints which should be imposed on the axion mass m a and its 
coupling constants for fulfilling the existing experimental and astrophysical bounds. Antici- 
pating our results, we remark that the consistent parameter region corresponds to large val- 
ues of x = tan/3 and so sin/? ~ 1. Hence, the axion coupling constants to the electron 
is Qae = 5 x 10~ 10 /g~ , and for the axion-nucleon coupling constants we have the follow- 
ing estimates g ap = —1.6 x 10 _6 /<r" and g an = 0.6 x lO -6 /^" (see e.g. in ref. Q), where 
h = (/a/10 6 GeV). 

The conservative interpretation of the experimental constraints from the axion search |8| 
imply the lower limit on the axion-nucleon coupling constants g a j\r < 10~ 4 or so, which trans- 
lates as fe > 10~ 2 . For m a > 2m e , the somewhat stronger limit emerges from the reactor 
search of a — ► e + e~ decay. Namely, for m a = 1.4 — 2 MeV one has Jq > 0.1 [20]. The latter 
limit, however, is subject of various theoretical uncertainties, and in the following we take the 
most conservative lower limit f a > 10 4 GeV. The existing limits on the axion decay into two 



photons are also fulfilled in this case [21]. 

For the standard DFSZ axion, since it is very very light, the severe limits on the scale f a 
emerge from the stellar evolution. In particular, the astrophysical bounds on the white dwarf 
luminosity function, helium burning lifetime of the HB stars or helium ignition time in low 
mass red giants imply g ae < (2.5 — 5) x 10~ 13 , and hence f a > (1 — 2) ■ 10 9 GeV II. However, 



these limits do not apply to our axion which is quite heavy and so its production rate in the 
stellar cores, with typical temperatures T up to 10 keV, is suppressed by the exponential factor 
exp(— m a /T). For m a > 300 keV this factor becomes enough small to render our model safe 
as far as the stellar evolution limits are concerned. 

However, stringent bounds on f a emerge due to the observed neutrino signal from the SN 
1987A which indicates that the supernova core cooling rate due to axion emission at the typical 
time t = 1 s from the collapse should not exceed 10 51 erg/s or so |J. If the axion-nucleon 
couplings are enough large, g a p,g a n > 10~ 7 , the axions are trapped in the collapsing core and 
are emitted with the thermal spectrum from an axiosphere with a radius R ~ 10 km, and 
energy luminosity at t = 1 s can be estimated as L a ~ / 6 ' x 3 • 10 50 erg/s [22, 23|. Hence, 



for fe<3 the axion luminosity L a becomes smaller than 10 51 erg/s, and the total energy and 
duration of the SN 1987A neutrino burst should not be affected. This consideration is certainly 
applicable for an axion with mass of order MeV, and so we impose the upper limit f a < 3 ■ 10 6 
GeV or so.0 

In addition, the supernova bounds provide strong limits on the axion decay rates into the 
ordinary as well as mirror photons. For m a ~ 1 MeV these decays are rather fast and strongly 
restrict the allowed parameter space. 

As far as the ordinary matter is transparent for the mirror photons, the emission of the 
latter can lead to the unacceptably fast cooling of the supernova core. The energy lose rate 
L' related to the mirror photon emission from the core volume can be estimated as follows. 
For f a < 10 7 GeV the axions are strongly trapped inside a core of the radius R ~ 10 km and 
have a thermal distribution. The temperature T is a function of the radial distance r from the 
center, and it also slowly changes with the time t elapsed from the collapse. The decay rate for 
an axion with the energy E into mirror photons is {m a /E)T' , where r' = T(a — > "f'j') is given 
by eq. (p0|). Therefore, we have L' (t) = 4TTm a V f Q r 2 n a (r,t)dr, where n a = 1.2T 3 /-7r 2 is the 
axion number density. Taking a core temperature T of about 20 — 30 MeV, one can roughly 
estimate that L' y ~ T'm a (4:irR 3 /3)(l.2T 3 /it 2 ) ~ 10" 2 g%ra\. Hence, the condition L'^ < 10 51 
erg/s implies that the decay width T', for m a ~ 1 MeV, should not exceed few s _1 . In other 
terms, using eq. (|i~9|), we obtain roughly the bound m 2 / f a < 10 -7 MeV. 

By taking the typical temperature profiles in the supernova core at different time moments 
||, one can compute the value L'(t) = 47rr'm a f Q r 2 n a (r,t)dr more precisely. In addition, 
m a and g' can be expressed as functions of x, x' and f a according to eqs. ( |l4] ) and (18). In 



Fig. 2 we show isocontours (dash-dott) corresponding to L' = 10 51 erg/s at t = 1 s from the 
collapse, which in fact indicate a lower limit on the PQ scale f a for given x and x' . 

Another constraint on our model emerges due to the axion decay into ordinary photons. For 
f a > 10 4 GeV the flux of axions emitted from the collapsing core is quite large, and thus if the 
axions would decay into the visible photons in their travel to the earth, one would observe the 



gamma ray burst associated with the SN 1987A. However, the SMM satellite data [24] sets the 
stringent upper limit on the 7 flux from the SN 1987A. In particular, these data indicate that 
photon fluence for for the photon energy band E~ = 4.1 — 6.4 MeV integrated over first 10 s from 
the collapse, did not exceed one photon per cm 2 . Therefore, the only way to to suppress the 
gamma ray signal from is to assume that the dominant part of the axions have undergone the 
decay inside the star envelope, during the flight time t s = cR s ~ 100 s, where R s ~ 3 • 10 7 km 
is a radius of the SN 1987A progenitor star, the blue giant Sanduleak. Hence, for the relevant 
axion energies E ~ 2S 7 = 8 — 13 MeV, the exponential factor exp(— m a ^tott s /E) should be 
very small, where the total decay width r to t incorporates the decays into the ordinary and 
mirror photons (f2C|), and if m a > 2m e , also into electron-positron (|2~l|). Therefore, T^ ot should 

12 The region f a > 10 10 is also allowed by the SN 1987A M , however it is not of our interest in this paper. 



not exceed few seconds. 

More precisely, by taking that the axions emitted from the axiosphere have a thermal 
spectrum with a temperature which at the time moment t from the collapse estimated as 
T a — f% (t/10 ms)" 1 ' 4 x 6 MeV [23], one can compute the expected photon fiuence in the 



energy band E~ = 4.1 — 6.4 MeV and confront it to the SMM bound. The results for various 
values of the parameter x are shown in Fig. 2. Here the dash curves correspond to the SMM 
limit on the 7 fiuence and so they set an upper limit on f a for given values of x and x' . 

Thus, we see that there is a parameter space compatible with the existing experimental 
and astrophysical limits, which yields f a around 10 4 — 10 5 GeV and m a about 1 MeV. This 
parameter range can be accessible for the experimental testing at the future reactor or beam 
dump experiments. 

4 Planck scale corrections 

A potential problem for the invisible axion models arises from the Planck scale effects. Many 
arguments suggest that the non-perturbative quantum gravity effects can induce the higher 
order effective operators cutoff by the Planck scale Mp which explicitly violate the global 
symmetry. Let us take the DFSZ model H] in which C7(1)pq symmetry is spontaneoulsy 
broken by the gauge singlet scalar S, (S) = f a /\^2 and consider a dimension 5 operator 
(XS(S'S) 2 /Mp + h.c, with A being some complex coupling constant. This term explicitly 
breaks the PQ symmetry and induces the mass term to the axion field: 

Ml = -^- * A/| x 0.17 GeV 2 (22) 

which should be compared with the dynamical mass term to 2 (||). The potential for 9 = a/f a 
then becomes f~ 2 P(a) = to 2 , cos # + M 2 cos(# + 5), where 5 is related to the complex phase 
of the coupling constant A and it is generically order 1. Although the first "genuine" term in 
this potential would fix the axion VEV at (9) = 0, the second term does not care about the 
strong CP problem and tends to minimize the potential along (9) = —5. Therefore, 9 < 10~ 10 
demands that M a < 10~ 5 m a , which in turn translates into the upper limit f a < A -1 ' 5 x 10 
GeV. This is in sharp contradiction with the lower limit f a > 10 10 GeV unless the constant A is 
extremely small, A < 10 -45 . Therefore, if one takes the Planck scale induced effects seriously, 
there is no room left for the invisible axion models. 

As for our model, it appears to be stable against the Planck scale corrections by threefold 
reason: (i) the PQ scale is lower, f a = 10 4 — 10 5 GeV, (ii) the axion is heavier, m a ~ 1 MeV, 
and (iii) the minimal possible explicit Planck scale operators in the Higgs potential are of 
dimension 6 and so they are suppressed Mp: Let us consider, e.g. the terms: 

\3 fu'lJi^ 



{H x H 2 f (H[H> 2 

1 2 P + Ml 

They induce the axion explicit mass term: 



x ~i^ +x -i^ + Kc - ^ 



Ml = 9X{X '^' 1)/fl4 ~ Ax'/e 4 x 0.5 • keV 2 (24) 

and the constraint M a < 10 -5 m a , for a typical value m a ~ 1 MeV, implies f a < 1.2(Aa;') _1 ' 4 10 5 
GeV, which is consistent with the parameter range found in previous section. 
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The Planck scale effects could explicitly break the PQ symmetry also in the Yukawa like 
operators fl2|| , e.g. 

atJ ^f d ^ Hl + aHJ ^f d U H ^ W 

which would interfere with the contribution to the Yukawa terms (||) and could deviate 6 from 
zero. However, for f a < 10 6 GeV these corrections are also negligible and do not give problems. 



Let us now comment about the model of Rubakov [17] which was based on the grand unified 
picture SU(5) x SU(5)' . It was assumed that the mirror GUT symmetry 577(5)' breaks down 
to the mirror standard model gauge group SU(3)' x SU(2)' x U(l)' at lower scale than the 
ordinary SU(5) breaks to standard gauge group. In other words, the VEV of the Higgs 24-plet 
$' of SU(5)' is somewhat lower than the VEV of the 24-plet Higgs of SU(5), ($) > ($'). Then 
the difference between the infrared poles for the strong coupling constants A' and A can emerge 
obtained due to the RG evolution from the corresponding scales down in energy. The gauge 
constant of SU(5) runs faster then the one of SU(3), and thus it is possible to have the bigger 
QCD scale in mirror sector, A' ^> A. In this way one could obtain obtain a very massive axion 
(e.g. m a ~ 100 GeV or more). 

However, in this case the M-parity is broken at GUT scales and therefore there is no strong 
reason to expect that it will be valid in the low energy Yukawa sector. For example, the Planck 
scale operators of dimension 5 are dangerous: 

G51077 + G'WlO'JF + £; -5-51077 + k'-^-WlO'JF... (26) 

M P M P y ' 

where H is the 5-plet Higgs. Due to M-parity we have G = G' and k = k! but the effective 
Yukawa coupling G = G + k(Q)/Mp is not equal to the mirror effective Yukawa coupling 
G' = G + fc($')/Mp. This should cause the big phase difference beteween them and so that 
9' ^ 9. Therefore, the model yjj is not stable against the Planck scale corrections. 

A way to overcome this problem is to break M-parity only at lower energies. E.g., one 
can take SU(5) x SU(5)' model with ($) = ($'), and break M-parity only due to different 
electroweak scales in two sectors, so that G' = G' . Clealy, this is nothing but embedding of 
our model considered in this paper in the grand unified picture. 

5 Conclusions 

We have suggested a new possibility for solving the strong CP-problem. It assumes that in 
parallel to the observable particle world described by the standard model, there exists a mirror 
world of particles with the identical Lagrangian, and two sectors share the same Peccei-Quinn 
symmetry realized d la Weiberg-Wilczek model. We asume that the mirror symmetry between 
two worlds is spontaneously broken so that the electroweak scale in the mirror sector v' is 
considerably larger than the ordinary electroweak scale v = 247 GeV. This in turn implies the 
infrared scale A' of the mirror strong interactions has to be somewhat larger than the ordinary 
QCD scale A ~ 200 MeV. In this way, the axion mass and interaction constants are actually 
determined by mirror sector scales v' and A', while the 9 terms are simultaneously cancelled 
in both sectors due to mirror symmetry. 

This means that the axion state essentially emerges from the mirror Higgs doublets, up to 
small (~ v/v') admixture from the ordinary Higgses. Hence, our axion is a WW-type axion 
for the mirror sector while from the point of view of ordinary sector it is DFSZ-type invisible 
axion with the symmetry breaking scale f a ~ v' S> v. 
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We have discussed the experimental and astrophysical limits on such an axion and have 
found an interesting parameter window where the PQ scale is f a ~ 10 4 — 10 5 GeV and the 
axion mass is m a ~ 1 MeV, which is accessible for the axion search in the future reactor and 
beam dump experiments. Interestingly, for such a parameter range our model is stable against 
the Planck scale corrections as far as the strong CP-problem is concerned. 

In our model, which in fact is the simplest mirror extension of the standard model, we 
cannot approach parameter space which could be applicable for explanation of the Gamma 



Ray Bursts via axion emission from the collapsing compact objects [25]. This, however, could 
be easily achieved in the supersymmetric version of our model, due to stronger effects in the 
renormalization group evolution of the strong coupling constants in two sectors. 
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